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SUMMARY 


2 * 0 41 ^ 

For initial stagnation pressures from 1 to 1,000 atm and stagnation 
enthalpies from 400 to 20,000 Btu/lb, nozzle-flow properties for equilibrium 
carbon dioxide have been computed and plotted on charts. Properties charted 
as a function of Mach number are as follows: temperature, pressure, density, 

speed, area ratio, dynamic pressure, stagnation-point pressure coefficient, 
Reynolds number, isentropic exponent, and molecular weight ratio. Tempera- 
tures, pressures, and densities across normal shock waves are also charted, 
and weight-flow rate is plotted as a function of stagnation enthalpy. 


INTRODUCTION 



Recently there has been much interest in the problems associated with 
flight to and within the atmospheres of Mars and Venus and return to Earth. 
Since the atmospheres of these planets are believed to be mostly carbon 
dioxide and nitrogen (e.g., ref. l) , experiments in hypervelocity tunnels with 
these gases and with air may yield important comparative information. 

Although some recent estimates indicate that the carbon dioxide content in 
the Venusian atmosphere may be only 20 mole percent (ref. 2) or even as low 
as 6 percent or less (ref. 3), other estimates have indicated an atmosphere 
of almost all carbon dioxide (e.g., refs. 4 and 5) • There is also consider- 
able uncertainty concerning the Martian atmosphere. Hence, until the chemical 
compositions of the atmospheres of these planets are better established, pre- 
liminary tests with both carbon dioxide and nitrogen alone and in several 
combinations are being considered. In order to define the conditions for 
tests in hypervelocity nozzles with these gases, various flow properties must 
be determined. In the present investigation equilibrium flow properties for 
carbon dioxide in hypervelocity nozzles have been computed. 

It is the purpose of this report to present /pharts of the computed 
nozzle-flow properties of carbon dioxide similar to those presented for air 
in reference 6 . The charts for carbon dioxide are for a range of stagnation 
pressures from 1 to 1,000 atm and stagnation enthalpies from 400 to 20,000 
Btu/lb. Properties charted as a function of Mach number are as follows: 
temperature, pressure, density, speed, area ratio, dynamic pressure, 
stagnation-point pressure coefficient, Reynolds number, isentropic exponent, 
and molecular weight ratio. Temperatures, pressures, and densities across 
normal shock waves are also included. Weight-flow rate is plotted as a 
function of stagnation enthalpy. Because in sin actual nozzle the flow will 



probably not be close to equilibrium for high stagnation enthalpies and low 
stagnation pressures, care should be exercised in the application of these 
charts . 


NOTATION AND CONSTANTS 


nozzle cross-sectional area, ft 2 
isentropic speed of sound, ft/sec 
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enthalpy (h = 0 for molecular gas at 0° R) , Btu/lb 
enthalpy, dimensionless 
characteristic length, ft 
Mach number 

molecular weight of mixture, lb/mole 

molecular weight of undissociated gas, 44.011 lb/mole for C0 2 
pressure, atm unless specified 
dynamic pressure pu 2 ^j , lb/ft 2 

universal gas constant: I .987 cal/mole °K; 8.3144 j/mole °K; 

1.987 Btu/mole °R; 15*5 ; or 0.7302 

Reynolds number, 

22.20 Btu/lb for C0 2 

entropy, Btu/( initial mole)(°R) 

entropy, dimensionless 

temperature, °R or °K as specified 

reference temperature, 492° R = 273° K 

speed, ft /sec 

weight-flow rate, Ib/sec 

molecular weight ratio, n^/m 
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isentropic exponent , ( 1-iH £ 

\o Zn p„ 

coefficient of viscosity, lb-sec/ft 2 
density, slugs/ft 3 

reference density, 0.00384 slug/ft 3 for C0 2 
sonic point 


Subscripts 

condensation 

reference condition 

reservoir or total condition 

conditions in front of normal shock wave 

conditions behind normal shock wave 


Conversion 

From Units In This Report To SI Units 1 



To convert from 


To obtain 

Physical quantity 

report units 

Multiply by 

SI units 

Area 

ft 2 

9.290X10 -2 

m 2 

Density 

slugs/ft 3 

5-i54xl0 2 

kg/m 3 

Dynamic pressure 

lb/ft 2 

47.88 

N/m 2 

Enthalpy 

Btu/lb 

2 . 324x10 s 

j/kg 

Length 

ft 

0.3048 

m 

Speed 

ft/sec 

0.3048 

m/s 

Temperature 

°R 

5/9 

°K 

Viscosity coefficient 

lb-sec/ft 2 

47.88 

N-s/m 2 

Weight- flow rate 

lb/ sec 

0.4536 

kg/s 


THERMODYNAMIC PROPERTIES AND VAPOR PRESSURES 

The thermodynamic properties of carbon dioxide which were used for the 
present study were computed by H. E. Bailey of Ames Research Center. His 
results (ref. 7 ) are recorded on magnetic tape for use in nozzle and body 
flow-field calculations with an IBM "JO^/lOk-O direct couple system. The data 
cover a temperature range of 100° to 25,000° K, a density range of 10 -7 to 
10 3 times earth sea-level density, and a pressure range of about 10 -7 to 
10 5 atm. There is very close agreement between the Bailey data and those of 
references 8 and 9« The reader should be cautioned, however, that because of 

The "SI system" (international System of Units, NASA TT F-200) . 


3 




a difference in choice of enthalpy reference datum., the value of h/RT 0 = 173-1 
should be added to the enthalpy results in reference 9 to make them compatible 
with those of reference 8 and Bailey (ref. 7) • 

Because it is important to avoid condensation conditions when testing 
models in a nozzle flow, a graph of vapor pressure as a function of tempera- 
ture is very useful. For ready reference vapor pressures as a function of 
temperature for carbon dioxide are compared with those for air, oxygen, and 
nitrogen in figure 1. The vapor -pres sure data for carbon dioxide, oxygen, 
and nitrogen were obtained from references 8, 10, and 11. The vapor pressures 
for air were estimated by assuming a perfect liquid or solid solution of 
oxygen and nitrogen in equilibrium with a 20-percent oxygen, 80-percent 
nitrogen vapor . 


CALCULATION OF NOZZLE-TLOW AND SHOCK PROPERTIES 


The nozzle-flow properties as a function of Mach number were desired for 
use in wind-tunnel testing. For convenience in the computational procedure, 
however, the flow properties were computed as a function of the ratio of 
static to stagnation pressure, with the corresponding Mach number being 
computed. Stagnation conditions p^. and h^ at the reservoir were specified 
for each case considered. To compute the area ratio, Mach number, and 
Reynolds number, the density was determined from the thermal equation of state 

p = pmo/ZRT (l) 

For p in slugs/ft 3 , p in atm, and T in °R, 

p = 1.873 p/ZT (2) 

where 

p = p t (p/pt) 

An isentropic expansion of the gas from the reservoir to downstream stations 
in the nozzle was assumed, and the molecular weight ratio Z, the temperature 
T, and the enthalpy h = (h/RT 0 )RTo were determined from the thermodynamic 
properties. The speed was then computed from the energy relation 

u = 223-6 >J~ht - h (3) 

The area ratio was calculated from the continuity equation 

A/A* = p*u*/pu (4) 

where the weight-flow rate at the throat (w/A* = p*u*) was determined as the 
maximum value of pu from a plot of pu vs p/pt for given stagnation con- 
ditions pt and h-t . (An alternate method which gives nearly identical values 
of weight-flow rate through the throat is given in ref. 12.) The Mach number 
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was computed from 


M = — = 


u 


/ 7 z iL t 

s] mo 


For T in °R, the speed of sound a in ft/sec was given hy 


a = 


(<&) = 33.61 fym 
N / W S 


(5) 


(6) 


where the isentropic exponent 7 was determined from 


7 = (f-*S-E 
\5 Zn 


( 7 ) 


The Reynolds number parameter was calculated from 

Re _ up 
Lp t PP t 


(8) 


The effects of dissociation and ionization on p were disregarded, and values 
of p were computed by (see ref. 13) 


= 2.28X1CT 8 T 1/2 ( 1 




( 9 ) 


for T in °R. 

An iterative procedure was used to obtain equilibrium temperature ratios, 
pressure ratios, and density ratios across a normal shock wave. With the 
Mach number Mx in front of the shock specified, the density ratio, p^/pi^ 
was assumed, and the pressure behind the shock was calculated from the 
expression 

P2 = Pi + (Pl^i) (Ui)(l - ^ (10) 

which resulted from combining the momentum and continuity equations. For the 
pressure in atmospheres 

Pa = Pi + ^ (Piui)(ui)(l - UD 

The local enthalpy behind the shock, h 2 , was determined by the energy equation 

h 2 = hi + 0 u^ - (1 (12) 
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For h 2 in Btu/lb 


1 


h 2 = hi 


50, 060 


UF 



(13) 


With p 2 and h 2 known, T 2 and Zg were obtained from the thermodynamic 
properties. The density behind the normal shock, p 2 , was then computed from 
equation (l), and Pz/Px was compared to the assumed value. This process was 
repeated until the assumed and computed values of Pz/Px were in agreement. 
The stagnation pressure and temperature behind the normal shock were then 
determined from the thermodynamic properties with an isentropic compression 
assumed from p 2 to p. , corresponding to an enthalpy change from h 2 to 


h t 2 = h ti* 


r'f > 
^2 


For use in computing forces and moments for entry-type capsules by 
modified Newtonian theory, the stagnation-point pressure coefficient c p s tag 
was determined by combining equation (10) and the Bernoulli equation for 
incompressible flow to give 


^Pstag 
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RESULTS 


The procedure outlined in the previous section was used to compute nozzle 
properties for equilibrium flow for various initial stagnation pressures and 
enthalpies. The initial stagnation pressures for the calculations were 1, 10, 
100, and 1,000 atm, and the initial stagnation enthalpies ranged from ^00 to 
20,000 Btu/lb. All calculations were made on an IBM 709^/70^0 system with the 
thermodynamic properties recorded in tabular form on magnetic tape. The 
nozzle-flow properties determined were temperature, pressure, density, speed, 
area ratio, Mach number, dynamic pressure, stagnation-point pressure coeffi- 
cient, Reynolds number, isentropic exponent, molecular weight ratio, and 
weight-flow rate. Temperatures, pressures, and densities across normal shock 
waves were also determined. 

Except for weight-flow rate which was plotted as a function of stagnation 
enthalpy, all properties were plotted as a function of Mach number and are 
presented on charts which are indexed on page 13* In the calculation of the 
characteristics with the machine program input values of p^ and (s/R)^ were 
used. As a result the determined values of ht varied as much as ±50 Btu/lb 
from those desired. The specified values of ht on each chart have been 
rounded off to the nearest 100 Btu/lb. For convenience in studies of flows 
around bodies, the local isentropic exponent 7 was plotted as a function of 
enthalpy (chart ik) as well as Mach number (chart 13 ) • The saturated vapor 
line for C0 2 in figure 1 (p. 9 ) was used to obtain the plot of stagnation 
enthalpy at condensation ht c vs M (chart 17 ) . The curves on the other charts 
of this report have not been terminated according to the condensation levels 
indicated in chart 17 * because these saturated levels obviously represent the 
upper enthalpy limits for condensation. Testing at lower total enthalpies may 
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be permissible for certain applications. The reader should observe that 
symbols are used on most of the charts only to identify curves , not to indicate 
the many points computed for each curve. 

Note that the plot of w/ptA* vs ht in chart 1 6 is extremely useful for 
deter min ing the stagnation enthalpy ht for a nozzle flow, since measurements 
of w, p t , and A* can be readily obtained. In reference 12 similar plots are 
given for various gases. 
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Figure 1.- Continued 
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(a) p > 10“ 1 atm 

Figure 1.- Vapor pressures of carbon dioxide, oxygen, nitrogen, and air. 
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Figure 1.- Concluded. 
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Chart 1.- Variation of temperature with Mach number. 
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Chart 2.- Variation of total temperature behind a normal shock wave with Mach number 
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Chart 3 • " Concluded 






Chart 4.- Concluded 
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Chart 5«- Variation of total pressure ratio across a normal shock wave with Mach number. 
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Chart 5*- Continued. 
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Chart 5*- Concluded 
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Chart 6 .- Variation of density parameter with Mach number; p = O.OO384 slug/ft 3 . 
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Chart 6.- Continued. 
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Chart 6.- Concluded. 
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Chart 7.- Variation of density ratio across a normal shock with Mach number 
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Chart 7 • ~ Cont inued. 
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Chart 7.- Continued. 
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Chart 7 * _ Concluded 
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Chart 8.- Continued 
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Chart 8.- Concluded.. 
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Chart 9*“ Concluded. 
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Chart 10.- Continued. 
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Chart 10.- Concluded. 






Chart 11.- Variation of stagnat 
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Chart 11.- Concluded 
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Chart 12.- Variation of Reynolds number parameter with Mach number. 
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Chart 12.- Continued. 






















Approximate velocity , u x 10" 



Chart l4.- Isentropic exponent as a function of enthalpy for various pressures 







Chart 15*- Concluded. 
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Chart 1 6 .- Variation of weight-flow parameter with total enthalpy. 
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NASA- Langley, 1965 A~10l6 
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Chart 17.- Total enthalpy for condensation-free flow, based on saturated vapor pressures in figure 





